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We develop a functional renormalization group approach which describes the low-energy single- 
particle properties of the Anderson impurity model up to intermediate on-site interactions U < 15A, 
where A is the hybridization in the wide-band limit. Our method is based on a generalization of 
a method proposed by Schiitz, Bartosch and Kopietz [Phys. Rev. B 72, 035107 (2005)], using two 
independent Hubbard-Stratonovich fields associated with transverse and longitudinal spin fluctua- 
tions. Although we do not reproduce the exponentially small Kondo scale in the limit U — > 00, the 
spin fluctuations included in our approach remove the unphysical Stoner instability predicted by 
mean-field theory for U > ttA. We discuss different decoupling schemes and show that a decoupling 
which manifestly respects the spin-rotational invariance of the problem gives rise to the lowest quasi- 
particle weight. To obtain a closed flow equation for the fermionic self-energy we also propose a new 
truncation scheme of the functional renormalization group flow equations using Dyson-Schwinger 
equations to express bosonic vertex functions in terms of fermionic ones. 

PACS numbers: 72.15.Qm, 71.27.-|-a, 71.10.Pm 



I. INTRODUCTION 

The Anderson impurity model (AIM) is one of the most 
important model systems in condensed matter physicSfi 
The model was originally proposed by AndersonS to de- 
scribe the properties of local moments in metals. Its 
Hamiltonian describes a single correlated impurity which 
is coupled to a band of non-interacting conduction elec- 
trons, 

a 

+ J2(^kdlck. + VkciJ.). (1.1) 

ktT 

Here, annihilates a non-interacting conduction elec- 
tron with momentum fc, energy dispersion and spin 
projection a, while the operator d^ annihilates a localized 
correlated d-electron with atomic energy and on-site 
repulsion U. The hybridization between the d-electrons 
and the conduction electrons is characterized by the hy- 
bridization energy Vk- We have also included in Eq. (jl.ip 
the Zeemann energy h associated with an external mag- 
netic field. 

The thermodynamic and spectral properties of the 
AIM can be calculated accurately by means of Wilson's 
numerical renormalization group- '4'^ (NRG), see Ref. ^ 
for a recent review. However, the calculation of the spec- 
tral function A{uj) of the rf-clcctrons by means of the 
NRG requires some computational effort, in particular 
if one needs accurate results for arbitrary frequencies uj. 
Since the quantitative knowledge of the spectral function 
of the AIM for all ui is essential in the context of the dy- 



namical mean-field theory describing strong correlations 
in realistic three-dimensional fermion systems,— it is im- 
portant to develop approximate analytical methods for 
calculating the spectral function of the AIM. Although 
several analytical approaches have been proposed to de- 
scribe the strong coupling regime^'^ d^d^d^d-^d'^ a satis- 
factory analytical alternative which can compete with the 
NRG in the strong coupling regime has not been found. 
This has motivated us to develop a new functional renor- 
malization group (FRG) approach to the AIM which is 
not based on the weak coupling truncation used by other 
authors. ^^"^^ Our approach extends the collective field 
FRG developed in Refs. PlUIii, which is based on 
the partial bosonization of the two-body interaction us- 
ing suitable Hubbard-Stratonovich fields. For simplicity, 
we focus here on the local moment regime^ where the 
energy Ed of the d-level is located below the Fermi en- 
ergy, but its double occupancy is prohibited by a strong 
on-site interaction U. In the limit U —>■ 00 the low- 
energy properties of the AIM are in this regime iden- 
tical to those of the Kondo model describing only the 
spin degree of freedom of the impurity. It is then nat- 
ural to decouple the interaction of the AIM in terms of 
Hubbard-Stratonovich fields representing collective spin- 
fluctuations, which should be treated non-perturbatively 
to describe the strong-coupling regime. 

Since we are interested in the correlation functions 
of the c?-electrons, we simply integrate out the conduc- 
tion electrons using the coherent state functional integral. 
The ratio of the partition functions with and without in- 
teraction at constant chemical potential /i and inverse 
temperature /3 can then be written as 

JV[d,d]e~So[d,d] ' ^1-2) 
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where the Gaussian part is given by 



So[d,d] = / dTy2d„{T)[dr+Ed- fi-<jh]d,{T) 

+ dT dr' VJ„(r)A^(r-T'K(r') 
Jo Jo „ 

= - /" ^ [zc^ - Co - A^M] J^^d^,, (1.3) 

J LJ _ 



with 

^0 ^ Ed ~ H ~ crh, 
and the interaction is 

/3 



(1.4) 



Su[d,d]=U dTd^{T)d^{T)d^{T)d^{T). (1.5) 

Jo 

Here, f = 4 denotes summation over fermionic 
Matsubara frequencies. Later we shaU take the hmit 
/3 — > CO where = / l^- The Fourier transform of 
the Grassmann fields d^ (t) in frequency space is defined 
by 



(1.6) 



(1.7) 



(1.8) 



J u 

and the hybridization function is 



with Fourier components 



A-(zc.) = 5] 



ioj — Efc + /i + (t/i 



The + 1-dimensional quantum field theory defined in 
Eqs. (|1.3p - (|1.5p will be the starting point of our further 
calculations presented below. 

Let us briefly outline the rest of this work: In Sec. [IT] we 
decouple the interaction (jl.Sp in the spin-singlet particle- 
hole channel using a complex Hubbard-Stratonovich field 
which describes transverse spin-flip fluctuations. By 
treating the resulting field theory on the level of the 
Gaussian approximation, we obtain results for the dy- 
namic structure factor and the self-energy in the ladder 
approximation. Of course, only for small U can we ex- 
pect the ladder approximation to give reliable results. 
For larger couplings we find a ferromagnetic Stoner in- 
stability, consistent with the Hartree-Fock approximation 
but in contrast to well-established results. To go beyond 
the simple ladder approximation, we use in Sec. IlIII a 
collective field FRG approach for our Bose-Fermi theory 
and implement a simple truncation to close the set of 
flow equations. Our inclusion of transverse spin fluctua- 
tions based on the FRG removes the unphysical Stoner 
instability and gives results which qualitatively describe 



the correct physics. However, although the quasiparticle 
weight vanishes with increasing interaction as it should, it 
does not obey the well-known Kondo scaling. To improve 
on these results, we introduce in Sec. llVI a new FRG ap- 
proach with partial bosonization in both the transverse 
and the longitudinal channel. We discuss the ambiguities 
with the distribution of weight among the two channels 
and show that in the manifestly spin-rotationally invari- 
ant case we get the strongest suppression of the quasipar- 
ticle weight. Finally, in Sec.|V]we summarize our results 
and discuss some open questions. 



II. LADDER APPROXIMATION IN THE 
SPIN-SINGLET PARTICLE-HOLE CHANNEL 

A. Hartree-Fock approximation and Stoner 
instability 

If we treat the interaction (|1.5p within the self- 
consistent Hartree-Fock approximation, we obtain the 
renormalized excitation energy = -|- 5(f , with 



(2.1) 



where we have used the notation a = —a for the spin 
label. The self-consistent Hartree-Fock Green function is 

1 



(2.2) 



(2.3) 



icj - - A'^(iw) 
The corresponding self-energy can be written as 

crm], 

where n = n-^ + is the average occupation and m — 
n^—ni is the average magnetization of the impurity level. 
For a sufficiently strong interaction, Hartree-Fock the- 
ory predicts a finite value of m even for h 0. In the 
wide-band limit, where the hybridization function can be 
approximated by 



A'^(icj) = -lA'^sgno;, 



(2.4) 



the Hartree-Fock self-consistency equation for the mo- 
ment m can be evaluated analytically at zero tempera- 
ture. For h = (where A"" = A is independent of the 
spin projection) and in the particle-hole symmetric case 
(where — fi = —U/2 and n = 1) the result is 

m — J ^ (tGq (iw) = — arctan ^^^^ . (2.5) 
This equation has a ferromagnetic solution (to 7^ 0) if 



Mo 



U 
ttA 



> 1. 



(2.6) 



This ferromagnetic Stoner- instability is an artefact of the 
Hartree-Fock approximationji which should be removed 
once fluctuation corrections are properly taken into ac- 
count. 
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B. Hubbard-Stratonovich transformation in the 
spin-singlet particle-hole channel 



In the strong coupling limit we expect that transverse 
spin fluctuations should play an important role to re- 
move the Stoner instability.^ It is therefore natural to 
decouple the interaction in the spin-singlet particle-hole 
channel with the help of a bosonic Hubbard-Stratonovich 
field, describing transverse spin fluctuations. To this end, 
we use the antisymmetry of the Grassmann fields to write 
the integrand in the interaction functional (jl.Sp as 

Ud^d^d^di = -U{d^d_Oidid^) = -Us{t)s{t), (2.7) 

where we have defined the composite spin-fiip fields, 

-sir) = dT(T)di(r) , sir) = d^{T)d^{T). (2.8) 

Introducing the complex bosonic Hubbard-Stratonovich 
fields xi''') x(t) conjugate to s(r) and s(r), the ra- 
tio (jl.2p of the interacting and non-interacting partition 
functions can be written as a functional integral over a 
six-component superfield = [df , d| , , d| , x, x] , 



(2.9) 



The Gaussian part of the bare action is 

J UJ ^ J UJ 

(2.10) 

and the interaction part can be written as 

= / [saXa + SaXa]- / ^SS," Lad^^a . (2.11) 

where 

so^ dT eJ^\s{T) ^ [ d^^d^+a,T, (2.12) 



dre-'^^sir)^ / d^+o,^d^i. (2.13) 



Because [Gq (ioj)]^-'^ in the Gaussian part (|2.10p of the 
action is by definition the inverse Hartree-Fock Green 
function (|2.2p . the Hartree-Fock self-energy 5^'^ should be 
subtracted from the interaction in Eq. (|2.1ip as a coun- 
terterm. The boson-fermion interaction in Eq. (|2.1ip can 
be written as 



[scsX^ + s^xM 



{uj + u), LO, uj)d^+i:^-\duj[XQ 

(uj - Ld,uj,U})d^^cdidu\Xo , (2.14) 



with the bare spin-fiip vertices given by 



-^(didix) 




{UJ-L0,UJ,U!) = 1. (2.15) 



A graphical representation of these vertices is shown in 

Fig.m 



FIG. 1: (Color online) Graphical representation of the bare 
boson-fermion vertices of the action Si[$] in Eq. (|2.1ip . Solid 
arrows pointing into the vertices represent d^, while outgoing 
solid arrows represent do- . The small arrows indicate the spin 
projections. The wavy arrows represent the spin-flip field x 
(pointing into the black dots) and x (pointing out of the black 
dots). The direction of the arrows associated with the wavy 
lines is chosen such that the incoming arrow adds spin to the 
vertex. 



Gaussian propagator of transverse 
spin-fluctuations 



If we integrate in Eq. (|2.9p over the fermion fields, we 
obtain the effective action S'cff[XiX] of the spin- flip fields 
X, X- Expanding Ses[Xi x] to quadratic order in the fields 
(Gaussian approximation), we obtain 



'S'cff[x,x] 



(2.16) 



The inverse spin-flip propagator in ladder approximation 
(LA) is given by 

[F^A{^iD)]-^^U-'-n^itQ), (2.17) 

where Hj]-(zcZ') is the non-interacting dynamic spin-fiip 
susceptibility 



(2.18) 



In the wide-band limit, where A'^^iuj) is given by 
Eq. (|2.4p . the integration in Eq. (|2.18p can be performed 
analytically at zero temperature. In this work we only 
need nf{^(i[D) in the absence of an external magnetic field 
and spontaneous ferromagnetism, where = ^ and 
A'^ = A are independent of the spin projection cr, 

-e + m+Ar 



TTld-Kld-l -t-2A) 



In 



(2.19) 



If we assume in addition particle-hole symmetry so that 
£.=£.o + S£_=0, then Eq. ((^19)) further simphfies to 



ttAHq (iw) 



ln[l 



k-|(i + k-|/2) 
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0{\ef), (2.20) 



where e = ui/A. Within the LA the spectral den- 
sity S'l^(w) of transverse spin-fluctuations (the dynamic 
structure factor) is then given by 



l-UU^{iQ) Jo 



(2.21) 
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FIG. 2: (Color online) Graph of the transverse spin struc- 
ture factor Sla(w) in ladder approximation for h — 0, uq = 
U/{ttA) — 0.9 and particle-hole symmetry. The characteristic 
energy scale of transverse spin fluctuations is in this approx- 
imation given by aj» « A(l — «o). 



or, equivalently, 

lmni^{LU + tO) - sgnc^^i^Adc^l). (2.22) 

A graph of S-^j^{u!) for uq — [//(vrA) = 0.9 is shown in 
Fig. [2] In the regime < 1 — wq 1 the LA predicts a 
weU defined peak in the dynamic structure factor at the 
energy scale 



= A(l — -uo)- 



(2.23) 



The width of the peak is of the order of A. In fact, in this 
regime the low- frequency behavior of S^j^{lu) can easily 
be obtained analytically. For ^ A we find 



1 



- Uo\LO\ 



SO that the corresponding structure factor is 

1 UJ 



(2.24) 



(2.25) 



For Uq —> 1 the energy scale vanishes, indicating an 
instability towards spontaneous ferromagnetism, as sug- 
gested by the Hartree-Fock approximation discussed in 
Sec. Ill Al Because we know that the AIM does not ex- 
hibit spontaneous ferromagnetism for arbitrary values of 
[/, we expect that the LA is only accurate for uq ^ 1. 
Before employing more sophisticated FRG methods to in- 
clude fluctuation corrections which remove this unphys- 
ical Stoner instability, it is instructive to consider the 
d-electron self-energy within the LA, which we shall do 
in the following subsection. 



D. Fermionic self-energy in ladder approximation 

The LA in the particle-hole channel gives the following 
self-energy for the d-electrons. 



'S"{iuj)= / F-^j^{iu!)GQ{iuj — iauj). 



(2.26) 



FIG. 3: Diagrams contributing to the fermionic self-energy 
E(iti;) in particle-hole ladder approximation, see Eq. (|2.26|l . 
Solid arrows denote fermionic Hartree-Fock Green functions 
and wavy lines denote the bare interaction. The first-order 
(Hartree) diagram is subtracted in Eq. (|2.27|l . 
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FIG. 4: (Color online) Graph of the quasiparticle residue Z as 
a function of uo = U /(ttA) within LA for h = and particle- 
hole symmetric filling Ed — 11= — y . While in the calculation 
of the solid line (blue) the full polarization was used, in the 
calculation of the dashed line (black) we have used the lin- 
earized spin susceptibility -k ISH-^ {iu)) ?a 1 — |aj|/A. 



The infinite series of Feynman diagrams which is included 
in Eq. (|2.26p is shown in Fig. [31 Subtracting the Hartree- 
Fock correction 5^'^ given in Eq. (|2.ip which formally 
arises from the counterterm in Eq. (|2.1ip , the frequency- 
dependent part of the self-energy can also be written as 

SYs^iioj) = T^^iioj) - SC I Ul^Aii^)G^{iuj - icrCj). 

J UJ 

(2.27) 

Of particular interest is the quasiparticle residue Z'^, 
which is defined by 



1 



1 



duj 



(2.28) 



In Fig. [4] we show the predication of the LA for the 
interaction-dependence of the quasiparticle residue in 
the non- magnetic, particle- hole symmetric case with 
and without linearization of the spin-fiip susceptibility 
IlQ{ioj). Obviously, the linearization of the spin-fiip sus- 
ceptibility does lead to modest changes of the Z-factor. 
The vanishing of Z at mq = 1 is an unphysical artefact 
of the LA, which implicitly contains the Stoner instabil- 
ity via the spin susceptibility. In the following section we 
shall show how to remove this instability by including the 



feedback of the fermionic wave function renormalization 
on the spin susceptibihty. 



III. FRG WITH PARTIAL BOSONIZATION IN 
THE SPIN-SINGLET PARTICLE-HOLE 
CHANNEL 

To go beyond the ladder approximation, we now use 
the collective field FRG approach developed in Refs. (TtI . 
[TsUToj to study the mixed Bose-Fermi theory defined by 
the action S[^] = 5o[$] + S'i[*] in Eqs. (I^IB and dTTTl) . 
In fact, this theory has a formal similarity with the the- 
ory describing two coupled metallic chains, which was 
studied by means of the FRG in Ref. [l^l- The exact 
FRG flow equations in the present problem are therefore 
represented by the same Feynman diagrams as given in 
Ref. 



A. Exact FRG flow equations 
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FIG. 5: (Color online) Exact flow equation for the fermionic 
self-energy T.'^{lu) with Hubbard-Stratonovich decoupling in 
the spin-singlet particle-hole channel. In our cutoff scheme 
only the bosonic propagator is regularized via a cutoff. Thick 
black arrows denote the flowing fermion propagator given in 
Eq. (|3.8|l . while thick wavy arrows with a slash denote the 
flowing single-scale spin-flip propagator deflned in Eq. (|3.5|l . 
The irreducible vertices are denoted by shaded triangles with 
the appropriate number of external legs. The dots over the 
fermionic two-point vertices on the left-hand side represent 
the derivative with respect to the RG cutoff A. 



To derive exact FRG flow equations, we modify our 
original model by introducing an infrared cutoff A into 
the Gaussian part (|2.10p of the action. In our renormal- 
ization group (RG) scheme the cutoff is introduced only 
in the bosonic sector i^^i^^i^^i^^ In previous FRG studies 
of the AIM, fermionic degrees of freedom were directly 
integrated out and a sharp Matsubara frequency cutoff 
was employed fi^ii^ For our purpose it is better to work 
with a smooth cutoff in order to avoid artificial singulari- 
ties introduced by a sharp cutoff. Formally, we introduce 
the cutoff via the following substitution in the bosonic 
part of the Gaussian action 5*0 [$] given in Eq. (|2.10p , 



where 



u- 



RA(iLo) 



C/-l+i^^(^cD), 



A 



■RM/K), 



(3.1) 



(3.2) 



and the function R{x) is given by^i 

R{x) = {l~ x)Q{l- x). (3.3) 
The flowing spin-flip propagator is then 



IJ_ 

" l + U[RA{iuj)-I{i{iLo)\ ' 



(3.4) 



where 11^ (iw) is the flowing irreducible transverse spin 
susceptibility. Introducing the corresponding single-scale 
propagator 



(3.5) 



the exact FRG flow equation for the irreducible self- 
energy of spin-(T electrons can then be written as 



djyT.l{lL0) 



Fj^ {iuj)G'J^{iuj — iauj) 



(3.6) 



Here, F 



(d„d„xx) 
A 



{uj,LO]Lo,Lij) is the flowing irreducible ver- 



tex with two bosonic and two fermionic external legs, 
and for the labels of the flowing three-legged boson- 
fermion vertices we have used the short notation = 
X and Xi — X: so that the non-zero combinations 



are F 



{d-TdiX-i'. 
A 



(w, to — uj^uj) — r 



A 



{lu, lu — uj,u!) and 

{lu — LU, uj, uj). Actually, 
the three-legged boson-fermion vertices have the symme- 
try 

ri'^"^'^^^)(^ - u,co,Cu) = r^j^^'^^\u;,co~ (3.7) 



so that the two vertices in Eq. (|3.6[) have the same value. 
The flowing d-level Green function G'J^{iuj) is related to 
the flowing self-energy via the Dyson equation, 



1 



■ eo--A-(iw)-i]X(*c^)- 



(3.8) 



A graphical representation of the exact FRG flow equa- 
tion p.6p is shown in Fig. [51 The initial values of the 
boson-fermion vertices are equal to the bare vertices 
given in Eq. (|2.15p . which is unity with our normaliza- 
tion. In principle, the integration of the FRG flow equa- 
tion (j3.6p generates also the Hartree-Fock contribution 
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(12. 3p to the self-energy as a boundary term which is of 
first-order in U and appears when we integrate this equa- 
tion up to A = 0. In practice, it is better to drop this 
first-order term and include the Hartree-Fock self-energy 
from the beginning into the rf-clcctron propagator, which 
amounts to imposing the initial condition 



i:i^{iuj) = SC = ^[n-am], (3.9) 



so that the initial G%^^{iuj) is the Hartree-Fock Green 
function. 



B. Truncation of the FRG equations via 
Dyson- Schwinger equations 



To obtain a closed system of RG flow equations, we 
need additional RG equations for the four-legged boson- 



fermion vertex F 

fermion vertices r^^''^'^'' and r^^'^^'^^ as well as for the 
flowing irreducible spin- flip susceptibility 11^ (iu)) which 
determines the single-scale spin- flip propagator F^{iiD). 
The exact FRG flow equations for the three- and four- 
legged boson-fermion vertices have been written down 
diagrammatically in Ref. [20^]; the crucial point is that 
the right-hand sides of these flow equations vanish at the 
initial scale A = Aq because they depend on higher order 
vertices which are not contained in the bare action. It is 
therefore reasonable to ignore the RG flow of the three- 
and four-legged boson-fermion vertices, which amounts 
to truncating the flow equation (|3.6p by replacing these 
vertices by their initial values, 



T^(dad„XX)(, , , ,. ,-, -,\ 



[UJ,U! - UJ,U!) 



Ao 



(w, LJ] LJ, LO) = 0, 



(3.10a) 

= 1, 
(3.10b) 

= 1, 
(3.10c) 



To close our system of flow equations, we still need 
an additional equation for the flowing spin-flip suscep- 
tibility I{j^{iu)), which in our cutoff scheme involves the 
pure boson vertex with four external legsi^ Fortunately, 
we can avoid the explicit analysis of this equation by us- 
ing the Dyson-Schwinger equation for our mixed boson- 
fermion theory, which implies an exact skeleton equation, 
relating the flowing spin-flip susceptibility to the flowing 
fermionic Green function G']^{iu}) and the flowing three- 
legged boson-fermion vertices. Using the same method 
as in Appendix B of Ref. [l3|, we obtain the skeleton 



FIG. 6: (Color online) Graphical representation of the skele- 
ton equation (|3.11|) relating the flowing spin-flip susceptibil- 
ity Iij^{iuj) to the exact fermionic Green functions and the 
three-legged boson-fermion vertices. The two diagrams are 
equivalent. 



for the three-legged boson- equation 



= G{iiLu)Gi{iuj - iu!) 

xTf''^'^\uj-u,uj,Lo). (3.11) 

A graphical representation of this equation is shown in 
Fig. El Note that Eq. p.7p guarantees that the two lines 
in Eq. p. lip are indeed identical. 

For simplicity, let us focus now on the non-magnetic 
case, where all correlation functions are spin-independent 
and we may omit the spin labels. Given our approxima- 
tion (j3.10bp . p.lOcp . the skeleton equation (|3.1ip then 
reduces to 

nf(iw) = - / GK{iuj)GK{iu -iu)), (3.12) 
while the FRG flow equation (|3.6p simplifies to 



5ASA(iw)= / F^{iuj)GA(iuj -iuj). (3.13) 

J UJ 

Eqs. p.l2p and p.l3p form a closed system of integro- 
differential equations for the flowing self-energy Yi\{iuj) of 
the d-electrons. Recall that Eq. (|3.13p depends implicitly 
on the flowing spin-flip susceptibility via Eqs. (|3.4p and 
p.Sp . Anticipating that there is no spontaneous mag- 
netism, the initial condition (|3.9p for the fermionic self- 
energy at scale A = Aq is simply 



Un 
~2~' 



(3.14) 



In the particle- hole symmetric case where Ed^^i = — C//2 
and n — 1 this cancels precisely the energy = Ed — I-il in 
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the Hartree-Fock d-electron propagator Gf^g{iuj), which 
in the wide-band hmit is therefore given by 



1 



iAsgncj 



(3.15) 



We next show that the solution of Eqs. (|3.12p . p.l3p does 
not suffer from the Stoner instability and correctly pre- 
dicts Fermi liquid behavior-i for arbitrary U . We shall re- 
strict ourselves to the particle-hole symmetric case from 
now on. 



the single-scale propagator defined in (13. 5p can then be 
written as 



e(A- 



^17 + ^' " 



(3.23) 



Actually, to be consistent with our approximation ()3.16|) 
which is based on the expansion of the fermionic self- 
energy to linear order in frequency, we should also expand 
the bosonic self-energy J{j:^{iu}) to linear order in uo. Using 



we have 



C. Low-energy truncation 

The coupled system of integro-differential equations 
given by Eqs. ([23), ([23]), (EH)), (|XT21) can be solved 
numerically for arbitrary frequencies, but this is beyond 
the scope of this work. Here, we shall focus on the low- 
frequency range \u}\ A. In this regime, it is reasonable 
to replace the flowing d-electron propagator appearing 
on the right-hand sides of Eqs. (|3.12p . p.l3p by the low- 
frequency Fermi liquid form 



GK{iuj) = 



Z, 



iAisgnoj 



(3.16) 



where the flowing wave function renormalization factor 
is defined by 



Zi 



1 



1 - 



diiuj) 



and the renormalized flowing hybridization is 
A, = Z,A. 



(3.17) 



(3.18) 



We consider Zi and A; to be functions of the logarith- 
mic flow parameter I — — ln(A/Ao). The approximation 
(I3.16P is certainly not sufficient at high frequencies, so 
that we cannot recover in this way the high-energy Hub- 
bard peaks at strong coupling which are known to appear 
at an energy scale of order [//2. 

Given the approximation (j3.16p . the frequency inte- 
gration in Eq. (|3.12p is easily carried out analytically, 
resulting in 



72 In 1 



1 + ^ 
^ + 2A, 



with 



ln(l 



(3.20) 



x{l +x/2) 

Introducing the dimensionless coupling constants 

ZiU 



Ul 



Ziua 
A 



ttA 
A 

a;' 



zfu 

ttA, ' 



(3.21) 
(3.22) 



TrAUiiiCo) « - y + 0{u!^). (3.24) 

As compared with the non-interacting susceptibility 
(I2.20p , the leading constant term in Eq. (|3.24p is reduced 
by the wave function renormalization factor, while the 
linear term is not renormalized. In this approximation, 
our single-scale propagator (|3.23p simplifies to 



TT e{A-\u\) 



1 



2 ' 



(3.25) 



which by construction of the regulator function given in 
Eqs. (|3.2p and (|3.3p is constant for frequencies below the 
running cutoff A. To determine Zi, we calculate the flow- 
ing anomalous dimension 77;, which is directly related to 
the derivative of the self-energy with respect to the flow 
parameter, 

r]i = -di InZi = ZiA lim -^dA^Aii^)- (3.26) 

Substituting Eqs. (|3.13p and p.l6p . we obtain 
d f FiiiQ) 



n = Z?A lim — - — - 

' ' t^^o d{iuj) 



iLu ~ iuj + iAisgii{uj — uj) 

(3.27) 

If we naively interchange the order of integration and 
differentiation, we encounter an ambiguous expression of 
the form 



hm . 
1 



1 



\UJ ■ 



~ iuj + iA;sgn(a; - 
2AiS{Q) 
Aisgnoj]^ ^ [A/(2e(w) - 1) 



u,) 



(3.28) 



where we have written sgncZ; = 2Q{lu) — 1. As pointed 
out by Morris, one should interpret the product of the 
delta function 5{x) with any function f{Q(x)) of the step 
function as 



(3.29) 



S{x)f{e{x)) = 6{x) / dtfit) 
Jo 

Using this relation, Eq. p.28p reduces to 



lim 



d 



d{iuj) 
1 

M + AiY 



1 



iuj — iuj + iA;Sgn(a; 
25{uj) 



(3.30) 
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n| , 1 , 1 , 1 , 1 

-10 -5 5 10 

FIG. 7: Numerical solution of Eq. (|3.33[) for different values 
of the bare coupling uo = U / (ttA). Here, is the scale where 
the running cutoff A is equal to the hybridization A, i.e., 
Aoe"'* = A. 




A I 1 1 , 1 , 1 , 1 , 1 

^0 2 4 6 8 10 
uq = [//ttA 

FIG. 8: The solid line is the wave function renormalization 
factor Z = limi_,oo Zi as a function of the bare coupling uq 
obtained from the numerical solution of our truncated flow 
equation (|3.33p . For comparison, we show as the dashed line 
the result of the ladder approximation. 



SO that Eq. (I3.27P takes the form 



(3.31) 



With the single-scale propagator given in Eq. (|3.25p , the 
frequency integration is trivial and wc finally obtain from 
Eq. (|3.27p for the flowing anomalous dimension 

The right-hand side of this expression depends on Zi via 
the flowing couplings ui and A; defined in Eqs. p.2ip and 
(P:^ so that 



diZi 



(3.33) 



is an ordinary differential equation for the flowing wave 
function renormalization factor Zi , which is easily solved 
numerically. The result for the flowing Zi is shown in 
Fig. [71 Obviously, for I — + oo the wave function renormal- 
ization factor approaches a finite limit, Z = lim/^oo Zi, 
which we show in Fig. [5] as a function of the bare cou- 
pling Uo- In contrast to the ladder approximation dis- 
cussed in Sec. |lT]we now obtain a finite Z for all values 
of the interaction, so that the fluctuations included in 
our simple FRG approach are sufficient to remove the 
unphysical Stoner instability. On the other hand, quan- 
titatively our truncation of the exact FRG flow equations 
does not reproduce the correct strong-coupling behavior 
of the quasiparticle residue, which is known to exhibit the 
same exponential suppression as the Kondo temperature 
for Uq ^ oo (see Ref. [23j). 



8uq 2 

exp[— TT Uo/ 

TT 



(3.34) 



From the numerical solution of our flow equation given 
in Eqs. (I3.32p . ()3.33|1 we find asymptotically 



0.445 



Uo > 1. 



(3.35) 



The fact that this is much larger than the exact result 
p.34p indicates that our approach underestimates the 
strength of fluctuations. In the following section we shall 
therefore include the longitudinal spin fluctuation chan- 
nel and show that it indeed further suppresses the value 
of Z, although within our approximations we are un- 
able to recover the exponential suppression described by 
Eq. I^M- 



IV. FRG WITH PARTIAL BOSONIZATION OF 
TRANSVERSE AND LONGITUDINAL SPIN 
FLUCTUATIONS 

A. Multi-channel Hubbard-Stratonovich 
transformation 

Using the antisymmetry of the Grassmann-fields, the 
local Hubbard interaction can be written in infinitely 
many equivalent ways, such as 



where 



m^ir) - U^s{t)s{t), (4.1) 



u. 



(4.2) 



Here, ^^(t) = d„{T)da-{T) and the composite fields 
Mt) = E^cr'^<T(T), s(t) = J|(T)dj.(T), and s(t) = 
di(T)d-\(T) represent the longitudinal and transverse spin 
components, see also Eq. (j2.8p . Note that by construction 
of the decomposition (|4.ip we have included both longi- 
tudinal and transverse spin fluctuations, but no charge 
fluctuations. In Sec. [TIT] we have satisfied Eq. (|4.ip by set- 
ting t/ll = and [/-'- = U, which is the natural choice if 
one is interested in bosonizing the transverse spin fluctu- 
ations. Alternatively, one could set C/" = U and C/^ — 0, 
so that longitudinal spin fluctuations can be introduced 



9 



via a suitable Hubbard-Stratonovich field. The ambi- 
guity associated with Eqs. (|4.ip and (|4.2p has been dis- 
cussed for many decades in the literature . ^^i^^'^^i^^'^^i^^ 
Depending on the physical problem of interest, certain 
special choices of C/H and can be advantageous At 
this point, we simply leave the precise values of C/" and 
U'^ unspecified but assume that both are positive and 
satisfy Eq. (|4.2p . Later we shall show that in the strong 
coupling regime of the particle-hole symmetric AIM the 
optimal choice is = U/3 and t/-*- = 2U /3. In this case 
the decomposition (|4.ip is manifestly spin-rotational in- 
variant and can be written as 

Un,{r)n^ir)^~^isir))\ (4.3) 

where the composite vector field s(t) — S {T)(jd{T) rep- 
resents the spin vector. Here, d'^ir) = [df (r), (i|(T)] and 
f? is the usual matrix vector of Pauli matrices. 

Starting from the representation (|4.ip , we decouple the 
interaction in the longitudinal spin channel using a real 
Hubbard-Stratonovich field rj and in the transverse spin 
channel using the complex Hubbard-Stratonovich fields 
X and X introduced in Sec. Ill Bl The partition func- 
tion can then be written as in Eq. ()2.9p . where $ = 
d|, dj, 7y, X, x] is now a seven-component super- 
field. For the Gaussian part of the action we then obtain 
instead of Eq. (|2T0l) 



iU'-)-'x^X^ + ^ I {U^^y'v^^r,^, (4.4) 



and the interaction part can be written as 



- f Y^SCd^ad^^. (4.5) 



The counterterm 5^'^ in Eq. (|4.5p is subtracted to correct 
for the inclusion of the Hartree-Fock self-energy in the 
propagator G^iiuS) in Eq. (|4.4p . While the transverse 
spin components and are given by Eqs. (j2.12p and 
p.lSp . the longitudinal spin component is given by 



r dTe'^^m{T) = [ Y^ad^^d, 

Jo Juj „ 



(4.6) 



For simplicity, we shall focus here on the non-magnetic, 
particle-hole-symmetric case. Since in the strong cou- 
pling limit the low-energy physics is expected to be dom- 
inated by spin fluctuations, we set up the ERG by intro- 
ducing a cutoff only in the bosonic fields 7], x ^^.d X as- 
sociated with spin fluctuations. Particle-hole symmetry 
guarantees that the exact self-energy at zero frequency 
cancels the term — /i, so that we may approximate 
Gl{iuj) by Eq. (IXTCD . 

B. FRG flow equations 

We now generalize Eq. p.ip by introducing cutoffs into 
both types of bosonic propagators appearing in the Gaus- 
sian action (|4.4p . 



(4.7) 



where a -L labels the two types of spin fluctuations, 
and the cutoff functions are 



R«{tLu) 
R^{iLj) 



2A ^f\uj\ 



7rA2 V A / ' 



7rA2 V A 



(4.8a) 
(4.8b) 



with R{x) given in Eq. (|3.3p . From the general FRG 
flow equations for mixed boson-fcrmion theories given in 
Refs. fl7, 19] we then obtain the following exact flow equa- 
tion for the self-energy, 



J to J UJ 

+ ( F^{iuj)Gl{iuj - iau))T^l'''^'''^''\ijj,uj- au),au))T^l'''^''''''\uj - au},u;,aiD), (4.9) 

J Q 

I 

which generalizes Eq. (|3.6p and is shown graphically in tion Xl ~ X s-nd xi — X-, see Eq. p.6p . The longitudinal 
Fig. m The superscripts label the different types of ver- and transverse single-scale propagators are 
tices; for the boson-fermion vertices associated with the 
transverse spin field we have used again the short nota- 

Fl(iu) = {~dt.Rl{iuj)\{Fl{iu:)\\ (4.10) 
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+ 




FIG. 9: (Color online) Graphical representation of the exact 
FRG flow equation (|4.9|l for the self-energies E\ and Ej^. The 
longitudinal spin-field rj is represented by dashed lines; the 
other symbols are the same as in Fig. [S] 



with 



F^itco) = [{U'^y -Illizu>) + RlitQ)]-\ (4.11) 



where n^(jcj) are the corresponding irreducible spin sus- 
ceptibihties. As discussed in Sec. IIIIBl we use the skele- 
ton equations (which can be formally derived from the 
Dyson-Schwinger equations) to relate the flow of the sus- 
ceptibilities to the flow of the d-electron propagators 
and the three-legged boson-fermion vertex. The skele- 
ton equation for the transverse susceptibility n^(i(D) has 
already been discussed in Sec. IIIIBl see Eq. p.lip : it can 
also be written as 



(4.12) 

The corresponding skeleton equation for the longitudinal 
spin susceptibility is 

n|(iw) ^ - I '^crGl{iuj)Gl{iuj - icD) 

xr(f''^"'')(w,w-w,w). (4.13) 

To obtain a closed RG flow equation for the self-energy, 
we still need flow equations for the vertices with three 
and four external legs appearing in Eq. (|4.9p . As in 
Sec. IIII Bl we simply set all mixed boson-fermion ver- 
tices with four external legs equal to zero, because these 
vertices vanish at the initial scale. However, in contrast 
to the FRG with only transverse spin fluctuations dis- 
cussed in Sec. IIII Bl the renormalization of the three- 
legged boson-fermion vertices is now important. In the 



approximation where mixed boson-fermion vertices with 
four and more external legs are ignored, the FRG flow 
equation for the longitudinal three-legged boson-fermion 
vertex is 

Fl{iu;')Gliiuj + + iid')Gl{iuj + iCu') 

Xi )^ (tJ -1- , (jj ) 

+ / Fj^{~iauj')G^f^{it^ + iuj + iuj')G^f^{iijj + i(:j') 



xT^^^'^'^'^uj + uj,0J + uj + u)\ -auj') 
y.Tf''''''\u^ + uj + Q',u: + uj',ud) 

y,Tf'd„'''\uJ+0j',UJ,-C7U^'), (4.14) 

while the corresponding flow equation for the transverse 
vertex reads 

Fl{iu>')Gl{iuj + + iLu')Gl{iuj + iCu') 
^j.{dM (ui + ij^u; + ij + Lj', -cj') 



Xr\ "(LJ +LJ ,UJ,UJ ). 



(4.15) 



Graphical representations of Eqs. (I4.14p and (|4.15p are 
shown in Figs. [TOl and [TT] From in Eq. (|4.5p we see 

that the initial conditions for the three-legged vertices at 
scale A = Ao are 

Tf;'^-''\uj + 10,10, to) = a, (4.16a) 

Y(f''^^^'\co + Ld,Lo,aLd) = 1. (4.16b) 



C. Low-energy truncation 

To make further progress, we now neglect the 
frequency-dependence of the three-legged boson-fermion 
vertices, setting 

Tf''"''^\co+LO,Lo,Lo) « ajl (4.17) 



-^idadiiXir) 
A 



(4.18) 



Keeping in mind that we have also neglected all vertices 
involving more than three external legs, our exact FRG 
flow equation (|4.9p for the self-energy reduces to 

^Am^io) = (7l)2 / Fliiio)GliiLO - iQ) 

J UJ 

+ [lif I Ft{iuo)Gl{ioo-iaLo). (4.19) 
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FIG. 10: (Color online) Graphical representation of the FRG 
flow equation (|4.14p for the longitudinal boson-fermion ver- 
tices r^'^T'^T'') (yppgj- graph) and r^'^^'*'^''' (lower graph). 
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FIG. 11: (Color online) Graphical representation of the FRG 
flow equation (|4.15|l for the transverse boson-fermion vertices 



(upper graph) and F 



A 



(lower graph). 



Moreover, the approximate FRG flow equation (|4.14p for 
the longitudinal spin-fermion vertex reduces to 

dAll = ill)' f_Fl{^u)[G\{^Lo)f 

- liilif I Ft{iCu)[Gi{~iCo)f , (4.20) 



while the FRG equation (|4.15p for the transverse spin- 
fermion vertex becomes 



^aTa - -7a(7;)' / Fjl{^u,)G^^i^Co)GX{^Co) . (4.21) 



Within the same approximation, the skeleton equations 
(I4.12P and (|4.13p for the spin susceptibilities are 



Jl 



nlitu) - -7I / 5]GX(*^)GX(*c.-z<i), (4.22) 
nj^{iu>) = -7a / Gl{iuj)Gl{iuj-iaLu). (4.23) 



Eqs. (I4.19p - (|4.23|) form a closed system of intcgro- 
differential equations for the frequency-dependent self- 
energy T,'J^(iuj) and the frequency- independent parts 7]! 
and 7a of the three-legged boson-fermion vertices. In or- 
der to make progress analytically, let us now approximate 
the flowing d-electron Green function on the right-hand 
sides of the flow equation by its low-energy Fermi liquid 
form dSjH). Then the integrations in Eqs. (l422l) . (|4^ 
can be performed analytically. For simplicity, we focus 
again on the non-magnetic particle-hole symmetric case 
and omit the spin-labels. We then obtain 



ni(.^) = ^-,^f (H I , (4.25) 



where the function /(x) is given in Eq. (|3.20p . The 
bosonic single-scale propagators can then be written as 



Fim = - 



e{A~\uj\) 



2Zf 



.,11 ^ A, II J \A,) 



e(A-|cD|) 



(4.26) 



(4.27) 



Here, we have introduced again the notation A; = 
A/A/ = A/(Z;A), and the running interaction constants 

m!' and uf- are defined by 



In terms of the dimensionless bare couplings Wq and Uq 
the condition (14. 2p reads 



u|}/2 + = uq. 



(4.29) 



As discussed in Sec. IIII C| for consistency we expand the 
function f{\uj\/Ai) in Eqs. (|4.24p and (|4.25p to linear or- 
der in frequency as in Eq. (|3.24p . Then the frequency in- 
tegrations in Eqs. (|4.19p . (|4.20p . (|4.2ip can be performed 
analytically and we obtain the following system of RG 
flow equations for the three running couplings Zi , jf , 
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and 7;-^ 



where 



diZi 
Zi 

drl 

if 
dilt 
it 



+ A 



I ' - li . 



[l/uf- + Xi - 



1 ' 



4-L 



^r = (7r)'A,/(i/< + A(-7r,7r 



(4.30) 
(4.31) 
(4.32) 

(4.33) 

(4.34) 
(4.35) 




Here, we have introduced the dimensionless integral 



FIG. 12: (Color online) Typical flow of the wave function 
renormalization factor Zi and of the frequency-independent 
parts 7;^ and 7]' of the boson-fermion vertices obtained from 
the numerical solution of Eqs. (|4.30|l " (|4.32|l . The scale = 
ln(Ao / A) is the same as in Fig. [7] The curves are for = 1 
and Uq — 1.5, such that uq — Uq/2 + uq = 2. 
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(4.36) 



Of special interest is the manifestly spin-rotationally 
invariant decomposition (|4.3p . In this case 

— U/3, such that each of the longitudinal and 
two transverse channels has the same weight. In terms 
of the dimensionless bare couplings ujj and Uq defined 



in Eq. (|4.28p we therefore have 



2uo/3. It is 



now easy to see that for any I we have uf = uf- = 2ui/3 

with ui — Ziuo and 7; = 7J' = 7^-^. Our flow equations 
(|i;3Dl)-(|132D then reduce to 



diZi 



Zi 
dill 
11 



rA, 



where 



m 



Au 

2 



2[3/(2uO + A( -7,] 
Ai = jfXiI {3/{2ui) + A, - 7,, 7, - 1, A,) 



(4.37) 
(4.38) 

(4.39) 
(4.40) 



The factors of 3 appearing here are related to the fact 
that for this choice of parameters we have three equiva- 
lent channels. Clearly, our truncated flow equations re- 
spect the spin-rotation symmetry of the problem such 
that for this special decomposition we can expect the 
best results. 



D. Results 

The three-dimensional system of differential equations 
given in Eqs. (|4.30p - (|4.32p can easily be solved numer- 
ically. An example of a typical RG flow is shown in 
Fig. 1121 Obviously, the wave function renormalization 
factor decreases monotonically as the RG is iterated, 
while the vertex corrections 7J' and 7/" both increase. 
The strongest variations occur at the scale U where the 
effective cutoff Aqc"'* is equal to the hybridization A. 

Before discussing the behavior of Z = limj^cxj Zi in the 
strong coupling regime, let us fix the optimal choice of the 
relative weight of the bare couplings Jj" and U-^ in the 
Hubbard- Stratonovich decoupling (|4.ip . In Fig. [T51 we 
show the dependence of the wave function renormaliza- 
tion factor on the choice of u^/uq for fixed values of uq. It 
turns out that all curves have a pronounced plateau with 
a minimum precisely at Uq/uq — 1, which corresponds 
to the manifestly spin-rotationally invariant decoupling 
(|4.3p . If we fix the choice of the Hubbard-Stratonovich 
decoupling by demanding that first-order variations of 
the results around the optimal choice should vanish, then 
we are naturally led to the decoupling (|4.3p . For the rest 
of this section we shall therefore use the manifestly spin- 
rotationally invariant choice Uq = Uq = 2uo/3. 

In Fig. [13] we show our numerical results for Z as a 
function of the dimensionless bare coupling uq- For com- 
parison, we also show the prediction of the ladder approx- 
imation, as well as our ERG results with decoupling only 
in the transverse spin channel, and numerically accurate 
results obtained within Wilson's numerical renormaliza- 
tion group. Our ERG calculation with simultaneous de- 
coupling in both transverse and longitudinal spin chan- 
nels obviously yields much better results for the suppres- 
sion of Z in the strong coupling regime than the single- 
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FIG. 13; (Color online) Dependence of the quasiparticle 
residue Z on the choice of the Hubbard-Stratonovich decou- 
pling, parameterized by the ratio ujl/Mj;!", with Uq/2-|-U(^ — uq. 
All curves exhibit a local minimum precisely at Uq/uq = 1, 
corresponding to the manifestly spin-rotationally invariant 
decoupling Uq — Uq — 2uo/3. While this minimum is also 
the global minimum for uq — 1,2 and 5, it is only a local 
minimum for uq = 0.5. In fact, for uq — 0.5 the curve is 
almost flat and has a global minimum at Uq/uq = and a 
local maximum between the two minima. 
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FIG. 14: (Color online) Wave function renormalization fac- 
tor Z = limi_,cx3 Zi as a function of the bare coupling uo — 
U / (vrA) obtained from the numerical solution of Eqs. (l4.3Up - 
(|4.32p . The two-channel FRG results are for the mani- 
festly spin-rotationally invariant choice of the two Hubbard- 
Stratonovich decouplings where u\ = Uq = 2uo/3. For a 
comparison, we have also shown NRG results from Ref. [l^ 
(which we have extrapolated to values uo > 4.5.) 



channel FRG discussed in Sec. lllll In fact, on the scale of 
Fig.fTHour FRG results for Z seem to be reasonably close 
to the exact numerical results, which in the strong cou- 
pling regime can be approximated by the asymptotic for- 
mula (|3.34p . To investigate whether our FRG approach 
reproduces the known exponential suppression of it 
is useful to present the data in Fig. [14] by plotting 1/Z 
on a logarithmic scale, as shown in Fig. (TSl Note that 
on this scale the exponential suppression of Z at strong 
coupling corresponds to a straight line. Obviously, for 
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FIG. 15: (Color online) Redrawing of Fig. 1141 the inverse 
wave function renormalization factor 1/^ is now plotted on a 
logarithmic scale. 




FIG. 16: (Color online) Low energy behavior of the spectral 
density A(lo) of the d-electrons as defined in Eqs. (I4.42p for 
uo = (dashed line) and 2, 4, 10. 



U > 15A our two-channel FRG results begin to devi- 
ate significantly from the NRG results and definitely do 
not reproduce the known exponential suppression of the 
quasiparticle weight for U — > oo. On the other hand, 
the two-channel FRG is accurate up to [/ « 15 A. We 
suspect that the deviations from the NRG results for 
U < 15A are due to our linearization of the spin sus- 
ceptibility, which is expected to lead to a suppression of 
the quasiparticle weight, see Fig. U) 

Finally, we show in Fig. [16] the spectral density of the 
d-electrons, which is defined by 



1_ 



-IniGiuj + iO). 



(4.41) 



For simplicity, we approximate G(aj-f iO) by G{^=^{iijj 
u) + iO) and obtain 



■nAA{Lo) 



1 



Z(w)A 



(4.42) 



with Z{uj) — 2^/=in(Ao/w)- The qualitative behavior of 
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A(uj) is in agreement with known results^: the width 
of the central Kondo peak is proportional to the wave 
function renormalization factor Z and its height at zero 
frequency is pinned to the value l/(7rA) in the particle- 
hole symmetric case considered here, so that the overall 
spectral weight of the low-energy peak is of order Z. The 
spectral line-shape at energies larger than A (including 
the broadened Hubbard bands at energy scales ±U/2) 
cannot be described within our low-energy truncation. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have proposed a functional integral 
approach to the Anderson impurity model which is based 
on the mapping of the original fermionic problem onto 
a mixed Bose-Fermi theory where the spin fluctuations 
are represented by bosonic Hubbard-Stratonovich fields. 
Our approach can be used to approximately calculate the 
spectral properties including the quasiparticle weight up 
to couplings U < 15A. 

This work also contains several technical advances 
which will be useful for other FRG calculations. In 
particular, we have shown that in FRG calcula- 
tions for Fermi systems using the technique of partial 
bosonizatiou j^^'^^i^^'^^i'^^'^^ the skeleton equations in the 
bosonic sector (which follow from the Dyson-Schwinger 
equations) can be used to close the FRG flow equations in 
the fermionic sector. Moreover, we have shown how the 
ambiguities inherent in multi-field Hubbard-Stratonovich 
decouplings can be resolved if one demands minimal 
sensitivity with respect to small variations of a given 
choice of decoupling. For a symmetric Anderson impu- 
rity model, this naturally leads to the manifestly spin- 
rotationally symmetric decoupling (j4.3p at strong cou- 
pling. With our decoupling scheme, this symmetry is 



also respected by the truncated flow equations. 
Our approach still has two major shortcomings: 

(i) We have not been able to reproduce the exponen- 
tial dependence Z oc exp[— 7r^ito/8] of the wave function 
renormalization factor for ?7 — > oo. Although we have 
tried several modiflcations of our approach, the strong 
fluctuations responsible for an exponential suppression 
of Z for U oo are apparently not correctly described 
within our truncation of the exact FRG flow equations. 

(ii) In order to make progress without using complicated 
numerical methods, we have made approximations which 
are only accurate at low energies < A. However, 
a direct numerical solution of our FRG flow equations 
should give rise to a much better agreement with NRG 
data at small couplings U and might also correctly repro- 
duce the high-energy behavior of the spectral function. 
Therefore one should numerically solve the coupled sys- 
tem of integro-diffcrcntial equations given by Eqs. (|3.12p . 
([XT^ or Eqs. (|iTg)) - (|i?^ . without approximating the 
fermionic Green function on the right-hand side of these 
equations by the Fermi liquid form (|3.16p . The fact that 
with our simple approximations we get results which 
agree reasonably well with NRG results gives us hope 
that solving this coupled system of integro-differential 
equations could lead to very accurate results. 
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